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Abstract. We discuss generation of magnetic field from cosmological perturbations. We consider the evolution of three com- 
ponent plasma (electron, proton and photon) evaluating the collision term between elecrons and photons up to the second 
order. The collision term is shown to induce electric current, which then generate magnetic field. There are three contribu- 
tions, two of which can be evaluated from the first-order quantities, while the other one is fluid vorticity which is purely 
second order. We compute numerically the magnitudes of the former contributions and shows that the amplitude of the pro- 
duced magnetic field is about ~ 10 _19 G at lOkpc comoving scale at present. Compared to astrophysical and inflationary 
mechanisms for seed-field generation, our study suffers from much less ambiguities concerning unknown physics and/or 
processes. 
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1. Introduction 

There are convincing evidences that imply existence of sub- 
stantial magnetic fields in various astronomical objects. Not 
only galaxies but systems with even larger scales, such as 
cluster of galaxies and extra-cluster fields, have their own 
magnetic fields (for a review on cosmological magnetic 
fields, see e.g., [I)). Conventionally the magnetic fields in 
galaxies, and possibly in clusters of galaxies, are considered 
to have been amplified and maintained by dynamo mecha- 
nism. However, the dynamo mechanism needs the seed mag- 
netic field and does not explain the origin of the magnetic 
fields. 

There have been many attempts to generate the seed 
fields. One of the approaches to this problem is to gener- 
ate magnetic fields astrophysically, often involving the Bier- 
mann mechanism j2). This mechanism has been applied to 
various systems: large-scale structure formation |3 1, ionizing 
front 1 4 1, protogalaxies 1 5 1 and supernova remnant of the first 
stars |6|. These studies showed the possibilities of magnetic - 
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field generation with amplitudes of order 10 16 ~ 10 21 G, 
which would be enough for the required field ~ 10 _20 G. 

On the other hand, cosmological origins, which are often 
concerned with inflation, can produce magnetic field with co- 
herence lengths of much larger scales, which is typically the 
horizon scale or possibly super-horizon scale f71l8ll9l [Tolll II . 
Also they can often produce fields with a wide range of length 
scales, while astrophysical mechanisms can produce fields 
only with their characteristic scales. For a constraints on mag- 
netic field with cosmological scale, see 1 12 1. 

Here, we consider magnetic-field generation from cos- 
mological perturbations after inflation. Around and after the 
decoupling, coupling among charged particles and photons 
become so weak that electric current can be induced by the 
difference in motions of protons and electrons. This elec- 
tric current leads to generation of magnetic fields. It is well 
known that vorticity of plasma produces such electric cur- 
rent El El El El El El- However, because vorticity is 
not produced at the first order in cosmological perturbations, 
we must study the second order. We will consider equations 
of motion for protons, electrons and photons separately up 
to the second order, although equation of motion for photons 
does not appear explicitly. To study electric current appro- 
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Fig. 1. All sky map of cosmological microwave background 
anisotropy obtained by the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite (upper) and schematic picture of 
cosmological magnetic fields generated from density fluctu- 
ations discussed here. Red (blue) regions are hot (cold) spots 
with a range of temperature ~ 2.725 K ±200/i K. The mag- 
netic field vectors are shown together with the map. Strong 
magnetic fields are generated by currents where the gradient 
of density perturbation in photons is large. 



priately, we must treat the three components separately 1 17 1. 
Furthermore, we will evaluate the collision term between 
electrons and photons, which is dominant and essential for 
the magnetic-field generation, up to the second order. From 
the equations of motion for protons and electrons, we obtain 
a generalized Ohm's law, which, combined with the Maxwell 
equations, leads to an evolution equation for magnetic field. 
Our study is based on the cosmological perturbation theory, 
which is highly successful in the anisotropies of the cosmic 
microwave background, and suffers from much less ambigui- 
ties concerning unknown physics and/or processes compared 
to astrophysical and inflationary mechanisms for seed-field 
generation. For details of our study, see 1 19l l20l . 



2. Formulation 

Euler equations for proton fluid and electron fluid are given 
by 



) \^p\v + ^EMp;i/) — + , (1) 



where T^Y*. and T^Mpte) ^ e ner gy- momen tum ten- 
sor of proton (electron) fluid and electromagnetic field cou- 
pling to protons (electrons) current, respectively. Here fi,u = 
0, 1, 2, 3 and i — 1,2, 3. The projection of the divergence of 
the energy-momentum tensors are computed as 

= (p + p)u'V. <1 + (S\ + u^JP'*, (3) 
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where p, P and j M are the energy density, pressure and elec- 
tric current, respectively. The r.h.s. of Eq. Q and 10 repre- 
sent the collision terms. Cpe — —Cip is the collision 
term for the Coulomb scattering between protons and elec- 
trons. This term leads to the diffusion of magnetic field and 
can be neglected in the highly conducting medium in early 
universe. On the other hand, the collsion terms for the Thom- 
son scattering between protons (electrons) and photons are 
expressed as C^yy- Because the collision term for the pro- 
tons can be neglected compared to that for the electrons, dif- 
ference in velocities of protons and electrons will be induced 
which leads to electric current. This electric current becomes 
a source for magnetic field. 

Now we evaluate the collision term for the Thomson scat- 
tering: 7(pi)+ e ~('7i) — > 7(p^)+e~(<7 2 '), where the quantities 
in the parentheses denote the particle momenta. The collision 
integral C\ e [/ (Pi)] f° r this scattering is given as 1 19 1, 
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where ot is the cross section of the Thomson scattering. Here 
moments of the distribution functions are given by 

— %PfliPi) = Pi, (6) 
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(2tt) 3 ^ ' ' ' 6 r 
where IE? is photon anisotropic stress. It should be noted that 
the collision term Q was obtained non-perturbatively with 
respect to the cosmological perturbation. 

Altogether, the Euler equations for protons and electrons 
are written as 

_ » m ,c p'-n (10) 



env£F* = 0, 



4<jTPyn 



(4 



(ID 



3 EVALUATION OF GENERATED MAGNETIC FIELD 



where m p is the proton mass, and the pressure of proton and 
electron fluids are neglected. We also assumed the charge 
neutrality: n = n e ~ n p . Subtracting Eq. i ll Ob multiplied 
by m e from Eq. il 11 multiplied by m p , we obtain 
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where w' 1 and j M are the center-of-mass 4-velocity of the pro- 
ton and electron fluids and the net electric current, respec- 
tively, defined as 



m e u^ 



(13) 



Noting the Maxwell equations Fff," = j M , we see that the first 
term in the l.h.s. of Eq. (1 1 2i are suppressed, compared to the 
second term, by a factor |21 1 
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where c is the speed of light, L is a characteristic length of the 
system and oj p = ^/ Anne 2 /m e is the plasma frequency. Sec- 
ond order vector perturbations are contained in the covariant 
derivative of the electric current and were evaluated in 1 16 1. 
They obtained the current magnetic field of order 10 _29 G at 
IMpc scale. As we will see below, this is much smaller than 
a value we obtain in this paper, which justifies the neglection 
of the first term in the l.h.s. of Eq. d!2i . 

The third term in the l.h.s. of Eq. (II 2i is the Hall term 
which can also be neglected because the Coulomb coupling 
between protons and electrons is so tight that \u z \ ^s> \u p — 
u\\. Then we have a generalized Ohm's law: 
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Now we derive the evolution equation for magnetic field. 
This can be obtained from the Bianchi identities Fi^x] = 0, 
as 

= e« fc u"Fy M 

= u»B% - ^ k C u\ + e** k C jtk 
u° ■ 

— {u % ,B 3 — u 3 ■ H £(B 3 u l - B l u 3 ), (16) 

u u 

where e ,J- * is the Levi-Civita tensor and B l = e l3k Fjk/2 is 
comoving magnetic field. We can expand the photon energy 
density, fluid velosities and photon anisotropic stress as 



(0) (1) (2) 

p 7 =P 7 + P y H , uo = 1+ «o H , 



(17) 



where the superscripts (i) denote the order of expansion. Re- 
membering that B l is a small quantity, we see that most of 



the terms in Eq. d!6i . other than the first and third terms, can 
be neglected. Thus we obtain 
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where the dot denotes a derivative with respect to the cosmic 
time, and we used the fact that there is no vorticity in the 
linear order: e t3k u = 0. The contributions of the first 
two terms in Eq. Jl 8I > were first noticed in 1171 . From this 
expression, we see that magnetic field cannot be generated 
in the linear order. Here it should be noted that the velocity 
of electron fluid can be approximated to the center-of-mass 
velocity at this order, u l e ^u 1 . 

3. Evaluation of generated magnetic field 

As we saw in Eq. Jl 81 . there are three main contributions 
to the generation of magnetic fields (1) baryon-photon slip 
term, (2) vorticity difference term, and (3) anisotropic pres- 
sure term. These terms derive from the fact that electrons are 
pushed by photons through Compton scattering when there 
exist velocity differences between them, or anisotropic pres- 
sure from photons. Here we now derive the power spectrum 
of magnetic fields, and then perform a numerical calcula- 
tion to evaluate it. The power spectrum of magnetic fields 

S(k) is defined as (\B{k)\ 2 ^ = S(k), where k is the wave 

vector, B(k) is Fourier component of magnetic fields. S(k) 
represents the expected variance of magnetic fields B(k), 
from which the component of the field with characteristic 
scale A can be derived through B\ w -J k 3 S (k) / (2ir 2 ) with 
A = 27r/fc. The power spectrum of magnetic fields S(k) is 
defined by the expected variance of the Fourier component of 
magnetic fields B(k) as S(k) = ^|B(fc)| 2 ^, where k is the 
wave vector. The component of the field with characteristic 
scale A can then be derived through B\ « y/ k 3 S(k) / (2tt 2 ) 
with A = 27r/fc. We consider a standard cosmological model 
which consists of photons, baryons, cold dark matter, neutri- 
nos, and cosmological constant, fixing all the cosmological 
parameters to the standard values. The density perturbations 
of them were solved numerically in a range of scales from 
10 kpc up to 10 Gpc, and they were then integrated to ob- 
tain S(k). We found that the field strength of generated mag- 
netic fields at cosmological recombination can be as large 
as 10 -16 8 G at IMpc comoving scale, and it becomes even 
larger at smaller scales (10~ 12 8 G at 10 kpc) (Fig. 0. Af- 
ter cosmological recombination, no magnetic fields would be 
generated, since most of electrons were combined into hydro- 
gen atoms and Compton scattering was no longer efficient. 
This means that the fields have an amplitude of 10~ 22 ' 8 G 
at IMpc (10~ 18 ' 8 G at lOkpc) at present because magnetic 
fields decay adiabatically as the universe expands after their 
generation. It is large enough to seed the galactic magnetic 
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fields required by the dynamo mechanism, which is typically 
of the order of 1(T 20 — 30 G at around the 10 kpc scaleU] 

Generated magnetic fields monotonically increase to- 
wards smaller scales on a range of scales in our calcula- 
tion. We found that the field has the spectrum S(k) oc fc 4 at 
scales larger than ~ 10 2 5 Mpc, which corresponds to super- 
horizon scales at recombination, S(k) oc k° at intermediate 
scales (10 2 ' 5 Mpc < A < 10 15 Mpc), and S(k) oc k 1 at 
scales smaller than ~ 10 Mpc, where the contribution from 
the anisotropic stress of photons dominates. This means that 
the field strength B is proportional to k 2 at the scales smaller 
than IMpc. If the primordial power spectrum of density fluc- 
tuations is given by a simple power law as predicted by infla- 
tion, our result implies that magnetic fields with the strength 
B w 10~ 12 ' 8 G unavoidably arise on 100 pc comoving scale 
at z ~ 10. This value is interesting for the evolution of struc- 
tures in the high-redshift universe, since those magnetic fields 
would be strong enough for magneto-rotational instability in 
the accretion disks surrounding very first stars (Population III 
stars) to be triggered, and affect the transport of their angular 
momentum 1 27 1 . The transport of angular momentum plays 
an important role for mass accretion process onto protostars. 
Since typical mass scale of the Population III stars is a key 
1 27 1 for the early reionization and chemical evolution of the 
universe, cosmologically generated magnetic fields should be 
one of essential ingredients in the model of the structure for- 
mation in the high-redshift universe. 

Since magnetic field creation mainly occurs when the 
modes of density perturbations with the corresponding scale 
enter the cosmic horizon and become causally connected 
1 20 1, the magnetic fields should exist at small scales below 
~ 10 Mpc even where the Silk damping effect by photons' 
diffusion has swept away the density perturbations at the last 
scattering epoch. Thus, in principle, the detection of magnetic 
fields below ~ 10 Mpc scales calculated here would tell us 
about density perturbations in photons (and baryons) in the 
early universe even at scales smaller than the diffusion scale 
at recombination. In this sense, the magnetic field generated 
in this mechanism can be regarded as a fossil of density per- 
turbations in the early universe, whose signature in photons 
and baryons has been lost. This result provides consequently 
a possibility of probing observations on how density pertur- 
bations in photons had evolved and been swept away at these 
small scales where no one can, in principle, probe directly 
through photons. 
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